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Abstract 

With  the  Combined  Release  and  Radiation  Effects  Satellite 
(CRRES)  measuremmts,  we  have  observed  an  extremely 
variable  outer  zone  relativistic  electron  population  from  25 
July  1990  to  12  October  1991.  Up  to  nov',  tUs  population  has 
been  modeled  by  the  static  NASA  solar  minimum  and  maxi¬ 
mum  models.  To  address  the  inadequacies  of  using  a  static 
model  to  describe  this  highly  dynamic  environment,  we 
develop  a  quasi-static  model  of  the  outer  zone  electrons  based 
on  the  readily  available  geomagnetic  activity  index,  Ap.  It  is 
^own  that  certain  quantities  used  to  parameterize  the  electron 
belt  morphology  are  moderately  correlated  with  the  logarithm 
of  the  15-day  running  average  of  Ap  (Ap,j).  We  therefore 
separate  and  average,  as  a  function  of  Ap,j,  the  438  daily 
average  radiation  belt  profiles  (electron  flux  versus  L)  for  each 
of  9  energy  chaimels  (1-8  MeV).  The  result  is  a  set  of 
average  flux  profiles  ^iiich  are  keyed  to  geomagnetic  activity. 
This  quasi-static  model  provides  a  more  accurate  r^resenta- 
tion  of  the  dynamic  outer  zone  electron  environment  than 
could  be  expected  from  any  static  model. 

I.  INTRODUCTION 

The  outer  radiation  zone  of  the  magnetosphere  is  populated 
by  trapped  relativistic  electrons  lying  between  2.4  and  8  Earth 
radii  (R^)  that  are  only  weakly  confined  to  the  magnetic 
equatorial  plane  [1].  It  presents  a  serious  radiation  hazard  to 
sensitive  components  in  satellites  and  may  limit  the  use  of 
emerging  technologies  in  space.  Compared  to  energetic 
proton  and  electron  populations  in  the  iimer  radiation  belt 
(within  a  radius  of  1.8  R^),  the  outer  zone  electrons  exhibit 
large  intensity  fluctuations  (several  orders  of  magnitude  in  flux 
change)  over  short  time  periods  (days)  [2,  3].  The  dynamical 
processes  that  drive  these  variations  are  still  not  identified,  but 
many  indicators  of  magnetospheric  storm  activity  show  some 
correlation,  albeit  weak,  with  outer  zone  intensifications  [4, 5, 
6].  Most  of  these  studies  have  been  conducted  with  data 
gathered  on  geosynchronous  satellites  flying  at  6.6  R^,  a  fixed 
altitude  at  the  outer  edges  of  the  outer  zone. 

The  high  variability  of  the  outer  zone  electrons  presents  a 
major  modeling  problem,  both  for  the  modeler  and  for  the 
user.  The  NASA  outer  zone  models  (AE8MIN  and  AE8MAX 
for  solar  minimum  and  maximum,  respectively)  were  prepared 
in  the  same  way  as  the  iimer  belt  models,  namely  long-term 
averages  applicable  for  missions  6  months  or  longer  [7]. 
Thus,  by  their  very  nature  they  are  inadequate  for  short 
missions.  Moreover,  recoit  comparisons  of  dose  measured  on 
orbit  with  predictions  from  the  NASA  models  indicate  that 
discrepancies  of  up  to  an  order  of  magnitude  exist  even  for 
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long  space  flights,  and  can  be  in  either  direction  dqiending  on 
shielding  thickness  [8,  9]. 

A  new  opportunity  to  measure  and  model  outer  zone 
electrons  is  presented  by  the  Combined  Release  and  Radiaticm 
Effects  Satellite  (CRRES)  which  measured  near-Earth  particle 
populations  over  a  wide  energy  range  from  25  July  19%  to  12 
CX;tober  1991.  The  CRRES  orbit  was  a  geosynchronous- 
transfer  orbit  with  an  18°  inclination.  Its  perigee  was  350  km, 
and  apogee,  33,500  km.  With  a  period  of  about  10  hours,  the 
satellite  made  at  least  four  transits  though  the  outer  zone  a 
day. 

In  this  paper  we  report  an  effort  to  model  the  outer  zone 
electron  population  in  a  frishion  that  gives  a  first  order 
estimate  of  its  dynamics.  We  modd  the  outer  bdt  over  the 
electron  energy  range  from  0.8  to  8  MeV  using  measurements 
from  the  High  Enmgy  Electron  Fluxmeter  (HEEF)  onboard 
CRRES.  We  first  establish  that  there  is  a  moderate  correlation 
between  the  electron  fluxes  and  the  IS-day  running  average  of 
the  global  geomagnetic  activity  index,  Ap,  delayed  by  one  day 
(referred  to  as  Ap,5).  The  Ap  index  is  made  readily  available 
by  the  NOAA-USAF  Space  Environment  Services  Crater  [10]. 
We  then  use  Api5  to  construct  eight  average  models  of  the 
electron  flux  variation  with  distance.  The  time  history  of  the 
outer  zone  can  then  be  roughly  reconstructed  for  any  time  for 
which  Ap  exists  using  the  eight  models.  The  model  fluxes  and 
those  of  the  NASA  AE8MAX  nxxiel  are  compared  to  the 
measured  profiles  to  determine  their  relative  accuracies  and 
limitations  in  predicting  the  outer  zone  high  energy  electron 
population. 

n.  INSTRUMENTATION  AND  DATA  HANDUNG 

The  High  Energy  Electron  Fluxmeter  (HEEF)  measures 
0.8  to  8  MeV  electrons  in  ten  differential  raergy  channels 
every  0.5  seconds.  The  instrument  is  a  telescope,  in  design, 
having  an  acceptance  half-angle  of  10.5°.  It  is  mounted 
perpendicular  to  the  spin  axis  of  the  satellite,  whi(±  in  turn, 
points  always  in  the  solar  direction.  The  satellite  spin  rate  is 
2  rpm.  Both  passive  and  active  (triple  coincidence)  shielding 
are  used  to  insure  that  high  energy  protoiis  and  low  rawgv 
electrons  do  not  contaminate  the  measurement.  The  instru¬ 
ment  and  its  extensive  pre-flight  calibration  are  described  in 
detail  elsewhere  [11].  The  results  presented  here  are  based  on 
a  preliminary  set  of  HEEF’s  calibration  constants.  A  set  of 
correction  factors  resulting  from  in-flight  calibration  analysis 
and  further  laboratory  tests  with  HEEF’s  backup  unit  is 
forthcoming. 

The  time  period  represrated  in  this  study  extends  from  27 
July  1990  to  H  October  1991,  about  14  months.  The  daubase 
is  generated  from  30  second  (1  spin)  average  HEF.F  fliiT#** 
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Figure  1.  Measured  (top  panel)  and  modeled  (bottom  panel)  daily  averages  of  1 .55  MeV  electron  flux  as  a  function  of  L  (in 
Rg)  and  time  (in  days). 


The  average  fluxes  are  accumulated  in  L-shell  (L)  bins  of  0.2 
Rc  over  the  outer  belt  region  from  2.4  to  6.6  Rg,  independent 
of  magnetic  local  time  and  magnetic  field.  The  Mcllwain  L- 
parameter  [12]  is  used  instead  of  altitude  because  it  identifies 
individual  drift  shells  confining  the  electrons  in  their  rotation 
around  the  Earth,  eliminating  the  need  to  sort  the  data  by  local 
time  or  longitude.  The  Olson-Pfitzer  tilt  dependent  magnetic 
field  model  [13]  is  used  to  determine  L,  and  although  the 
uncertainty  in  L  increases  for  L  >  4.5  during  magnetically 
active  periods,  the  model  provides  a  reasonable  basis  for 
ordering  the  data.  The  data  in  each  L  bin  is  averaged  to 
obtain  a  daily  flux  value  for  each  altitude  (L)  and  for  each 
energy  channel.  These  daily  spin-averaged  flux  values  are 
used  to  generate  the  models.  We  treat  the  spin  averaged 
fluxes  as  unidirectional  fluxes  with  units  of  e  (cm’ s  sr  keV) 
For  this  to  be  strictly  valid,  the  electron  population  must  be 
isotropic.  This  was  determined  to  be  approximately  true  after 
examining  a  database  which  is  a  superset  (includes  pitch  angle 
information,  and  has  a  higher  spatial  resolution^  of  the  one 
u.sed  in  this  study.  We  quantified  the  shape  of  each  pitch 
angle  distribution  by  defining  a  pitch  angle  index  as  follows. 
For  a  given  pitch  angle  distribution,  the  electron  flux  was 


averaged  over  both  65-90^  (j*,)  and  40-65"  ([55)  pitch  angle 
intervals.  A  pitch  angle  index  was  then  formed  from  the 
logarithm  of  the  ratio  j,o  ^  Js.'s  •  Thus,  a  perfectly  isotropic 
distribution  would  have  an  index  of  0.  An  automated  survey 
over  the  entire  1.5  MeV  electron  data  set  (restricted  to  L 
between  2,5  and  6.5)  revealed  that  81  %  of  pitch  angle 
distributions  had  an  index  between  -.3  and  +.3.  In  no  cases 
did  the  index  exceed  an  absolute  value  of  1. 

III.  Parameters  for  describing  and 
Ordering  The  Outer  Zone  Electron 
Flux  Variations 

The  measured  1.55  MeV  daily  average  electron  fluxes  over 
the  CRRES  mission  lifetime  are  shown  in  the  top  panel  of 
Figure  1  (the  bottom  panel  will  be  discussed  later).  The  flux 
values  (encoded  by  gray  scale)  range  from  <!()■'  (a  >  1(V  e 
(cm’ s  sr  keV) ',  and  are  plotted  as  a  function  of  L  (in  R^)  on 
the  vertical  axis  and  time  (in  days)  on  the  horizontal  axis.  The 
time  interval  between  tick  marks  is  10  days.  The  wide 
variation  in  the  radial  position  of  the  inner  edge  of  high  flux 
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regicms  and  the  overall  flux  levels  themselves  is  quite  evident. 
For  periods  of  more  intoise  fluxes,  the  inner  edge  is  located 
between  2.5-3  Rg  and  the  fluxes  peak  above  10*  e  (cm*  s  sr 
keV)  '.  By  contrast,  days  of  less  intoise  fluxes  have  an  inner 
edge  of  maximum  intensity  around  4.5  Rg  and  a  flux  intensity 
peak  below  10*  e  (cm*  s  sr  keV)  '.  To  illustrate  these  differ¬ 
ences  in  more  detail  we  show,  in  Figure  2,  flux  profiles  as  a 
function  of  L  for  8  August  1990  and  8  January  1^1,  times  of 
intense  and  weak  outer  zone  electrons,  respectively. 


We  call  Pj,  the  'centroid  L*  parameter,  defined  in  analogy  to 
the  colter  of  mass  as: 


P,  (E)  ' 


^  (£)  .Ai, 

i»l _ 

21 

(£■)  .Ai, 


Figure  2.  Eleciron  flux  inlennty  u  «  Amction  L*vitue,  in  Earth  ndii,  for 
an  intenae  outer  zone  (8  August  1990,  solid  line)  and  a  weak  outer  zone  (8 
January  1991,  dashed  line).  Arrows  mark  the  centroid  L  of  each  profile  (see 
text). 

N(rt  only  do  the  flux  peaks  differ  in  intensity  and  position,  but 
the  form  of  the  profile  is  extremely  different  in  the  two  cases. 
The  August  profile  rises  steqply  with  increasing  L  and  then 
falls  off  relatively  slowly.  The  January  profile  both  rises  and 
falls  off  slowly,  with  similar  slope  magnitudes.  To  model  the 
outer  zone  electrons  we  need  to  find  an  ordering  parameter  for 
these  different  profiles  of  electron  flux  intensity  as  a  function 
of  L.  To  do  fliis  we  first  characterize  the  profiles  with  two 
parameters  likoied,  for  a  mass  distribution,  to  the  total  mass 
and  the  center  of  mass. 

For  an  individual  flux  profile,  such  as  either  of  those 
shown  in  Figure  2,  for  a  given  energy  channel,  we  define  two 
parameters  P„  indicating  flux  intensity  in  the  outer  zone,  and 
P 2*  indicating  the  coiter  of  the  flux  distribution.  We  call  P,, 
the  'profile  flux*  parameter.  It  is  the  integral  of  the  flux  over 
L: 


21 
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In  these  equations: 

Lf  =  L-shell  bin  midpoint  in  R^  (2.S  to  6.5  R^); 
delta  L  =  L  bin  widA  (0.2  R£); 
i  =  L  bin  number  (1  to  21);  and 
ji(E)  =  flux  for  a  given  oiergy  duume!  (E)  and  L  bin  (i). 
In  Figure  2,  arrows  mark  the  crmtroid  L  for  die  two  profiles 
at  3.8  and  4.6  R^  for  the  intense  and  weak  profiles,  respective¬ 
ly.  The  corresponding  profile  fluxes,  P,,  are  1.4x10^  and  7. 3x 
10'  having  units  of  flux  times  R^. 


Figure  3.  Profile  flux  (P,;  lop  piiteO  tod  ceniroid  L  (P,;  middk  ptnel),  m 
defined  in  die  text,  plotted  venui  time  (m  dayi).  Both  panmeten  were 
derived  fi»m  daily  1 .55  MeV  fluxet.  Fifteen  day  running  avenge  of  Ap 
(A}>,])  vemii  time  (in  day*)  it  plotted  in  bottom  panel. 

The  two  parameters  defined  above  are  calculated  for 
electron  flux  profiles  for  each  day  of  the  CRRES  mission  and 
for  each  HEEF  energy  chaimel.  Figure  3  shows  these  two 
parameters  for  1.55  MeV  electrons  (P,,  top  panel;  Pj,  middle 
panel)  plotted  versus  day.  The  bottom  panel  is  a  15  day 
running  average  of  the  Ap  magnetic  activity  mdex  (referred  to 
as  Api,  for  the  remainder  of  this  paper).  Ap  is  a  readily 
available  index  of  magnetic  activity,  constructed  from  ground- 
based  magnetic  stations  having  a  broad  distiibution  in  latitude 
and  longitude  [10].  Ap  varies  linearly  with  the  magnitude  of 
the  magnetic  disturbance.  Figure  3  shows  that  magnetic 
activity  was  most  quiet  (Ap,,  ^  - 10  nT)  from  November 
1990  through  February  1991.  In  this  quiet  period  P,  (profile 
flux)  has  relatively  low  values  and  Pj  (centroid  L)  has  relative¬ 
ly  high  values.  In  general,  as  magnetic  activity  diminishes,  P, 
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(profile  flux)  decreases  and  (centroid  L)  increases.  During 
periods  of  high  magnetic  activity  the  reverse  is  true,  as  can  be 
seen  in  the  periods  before  December  1990  and  after  March 
1991.  These  trends  are  quantified  in  Figure  4  where  P,  (top 
panel)  and  Pj  (bottom  panel)  are  plotted  versus  log(Ap,j)  with 
an  imposed  1  day  lag  time.  A  lag  time  of  1  day  was  chosen 
because  it  was  found  to  minimize  the  error  between  the 
measured  and  modeled  fluxes  to  be  discussed  in  the  next 
section  of  this  paper.  A  linear  regression  of  the  437  data 
point  yields  a  correlation  coefficient  of  -  .65  for  both  the 
direct  relation  between  P,  and  log(Ap,j)  and  the  inverse 
relation  between  Pj  and  log(Ap,j).  The  eight  large  dots  on 
each  plot  refer  to  modeled  parameters  and  will  be  addressed 
later.  Given  the  moderate  correlations  between  the  flux 
profile  parameters  and  log(Ap,]),  we  use  Ap,j  to  sqrarate  and 
bin  the  electron  data  to  produce  our  flux  models. 


LOGC Ap  ts) 

Figure  4.  Profile  flux  (P,;  top  panel)  and  centroid  L  (Pj;  middle  panel) 
plotted  versus  log(Ap,j).  A  linear  regression  was  performed  on  the  437  data 
points,  with  the  best-fit  line  plotted  and  the  correlation  coefficient  given  in 
upper  right  corner  of  panel.  Eight  large  dots  are  the  profile  parameters 
derived  from  the  eight  model  profiles. 

IV.  Outer  Zone  Electron  Models 

Because  of  the  highly  dynamic  nature  of  the  outer  belt 
electrons,  a  single  model  to  reproduce  all  conditions  was 
deemed  impractical.  After  much  trial  and  error,  we  concluded 
that  8  models  divided  over  the  CRRES  Ap,,  range  between  5 
and  SS  nT  would  give  a  good  representation  of  most  condi¬ 
tions  measured  on  CRRES.  All  daily  flux  profiles  corre¬ 
sponding  to  the  days  assigned  to  one  of  the  8  Ap,,  intervals 
(assuming  a  1  day  lag  time)  were  averaged  together  to 
determine  the  model  profile  for  that  activity  level.  For 
example,  model  profile  #1  is  the  average  of  the  individual  flux 
profiles  for  the  31  days  in  which  the  value  of  Ap„  on  the  day 
preceding  the  daily  flux  profile  was  between  5  and  7.5  nT.  A 
sununary  of  the  model  inputs  is  given  in  Table  I.  Included  in 
the  table  is  the  range  and  the  corresponding  average  of  Ap,„ 
the  number  of  daily  profiles  used  to  compute  the  average,  and 
the  percent  of  the  total  days  (438)  u.sed  for  the  profile. 


Table  I  -  Ap,,  Model  Statistics 


Model# 

APa  <sT) 

Aw  Ap„ 

Dtfttfcn 

lUace 

(bT) 

Mcdd 

Toul 

1 

5.t^7J 

0.8 

3) 

7.1 

3 

7  10.0 

8.7 

S3 

lt.9 

3 

10.0-13.5 

11,2 

35 

8.0 

4 

12.5.13JJ 

13.7 

35 

8,0 

S 

lS.0-30.0 

17.J 

98 

22.4 

6 

30.0-2S.0 

22.4 

71 

16.2 

7 

2S.0-1S.0 

21.9 

56 

12J 

% 

35.0.55.0 

42J 

29 

6.6 

The  resulting  model  profiles  for  four  of  the  energy 
channels  are  shown  in  Figure  5  with  electron  flux  plotted 
against  L  shell.  For  the  sake  of  clarity,  only  every  other 
model  profile  is  displayed.  Model  profiles  #\  and  Wl  (#1  is 
the  bottom  curve  displayed  in  each  plot)  which  were  averaged 
from  the  extended  quiet  interval  from  November  1990  - 
February  1991,  are  significantly  different  in  shape  and 
magnitude  from  the  remaining  higher  activity  profiles.  The 
higher  activity  profiles  are  themselves  very  similar  in  both 
shape  and  magnitude.  To  show  how  well  the  models  can 
reproduce  the  original  data,  the  P,  and  P2  model  parameters 
from  the  1.55  MeV  chaimel  were  calculated  and  plotted  as 
solid  circles  on  their  respective  plots  of  Figure  4.  As  expect¬ 
ed,  the  model  parameters  P,  and  Pj  vary  with  log(Ap„)  in  the 
same  way  as  do  the  daily  profile  parameters.  This  relationship 
holds  true  for  all  energy  channels. 

The  model  flux  profiles  ordered  by  geomagnetic  activity 
shown  in  Figure  5  are  also  compared  with  the  single  flux 
profile  for  each  of  two  static  models  briefly  described  next. 
To  compare  the  new  Ap  dependent  models  with  existing  static 
models  and  data  sets,  the  NASA  solar  maximum  electron 
radiation  belt  model  (AE8MAX)  was  used  to  construct  a  flux 
versus  L  profile  for  each  energy  channel  of  Figure  5. 
AE8MAX  is  a  matrix  of  omnidirectional  integral  fluxes  stored 
as  a  function  of  threshold  energy,  B/Bq,  and  L  parameter  [7], 
B/Bq  is  the  ratio  of  the  magnitude  of  the  magnetic  field  at  the 
point  in  question  to  that  on  the  magnetic  equator  following  the 
same  field  line,  and  is  a  measure  of  magnetic  latitude.  It, 
therefore,  requires  some  preprocessing  for  a  direct  conqiarison 
with  CRRES  flux  data.  To  use  AE8MAX  to  predict  the 
electron  fluxes  CRRES  should  see,  CRRES  orbits  were  traced 
through  this  matrix,  and  the  omnidirectional  integral  fluxes 
converted  to  unidirectional  differential  fluxes  for  each  of  the 
CRRES  electron  energy  channels.  The  comparison  is  shown 
in  Figure  5  as  the  solid  curve.  In  all  energy  channels  the 
NASA  model  gives  values  higher  than  the  CRRES  models 
above  an  L  of  5  Rg.  For  energies  above  about  2  MeV  the 
NASA  model  gives  higher  values  at  all  L  values  above  —3.4 
Re- 
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Figure  5.  Eight  Ap„  model  profiles  for  4  different  energy  channels  (as  specified  in  key).  Also  included  are  the  NASA 
AE8MAX  model  (heavy  solid  line)  and  the  CRKES-static  model  (dotted  lines)  profiles. 


To  show  that  the  NASA  mode!  electron  fluxes  are  indeed 
higher  than  aU  the  CRRES  data,  and  to  illustrate  the  modeling 
benefit  of  quasi-static  models  relative  to  fully  static  models, 
we  also  cmted  a  single  static  model  from  the  CRRES  data. 
The  CRRES-static  model  is  sinqily  an  ordering  of  the  full  set 
of  HEEF  data  by  L  shell  (0.2  bins)  to  yield  one  average 
flux  profile  for  each  energy  channel.  The  electron  flux 
profiles  for  the  CRRES-static  model  are  also  included  in 
Figure  S  and  are  drawn  with  dotted  lines.  The  CRRES  static 
model  fluxes  are,  not  surprisingly,  similar  to  the  quasi-static 
model  fluxes  for  active  magnetic  conditions,  since  in  linear 
averages  of  quantities  varying  by  orders  of  magnitude,  the 
high  values  will  dominate.  The  profile  intmsity  and  shape  of 
the  average  CRRES  model  fluxes  are,  however,  still  very 
different  from  and  well  below  the  NASA  model  levels. 

V.  DISCUSSION  And  Conclusions 

Since  daily  Ap  magnetic  indices  are  disseminated  weekly 
[10],  the  quasi-static  model  presorted  here  provides  a  noethod 
for  calculating  the  short  term  average  fluxes  of  the  outer  zone 
electron  belts.  To  see  how  well  this  would  have  worked  for 
CRRES  while  in  flight,  we  used  the  daily  values  of  Ap,}  from 
Figure  3  to  identify  which  model  electron  profile  to  use  for 


each  day  of  the  CRRES  mission.  The  survey  plot  derived 
from  the  noodels  is  plotted  in  the  bottom  panel  of  Figure  1  for 
direct  visual  conqrarison  with  the  measured  values.  The 
agreem^t,  although  not  perfect,  is  far  better  than  any  static 
model  could  give.  The  models  maintain  the  overall  gross 
features  of  the  measured  data  such  as  the  quiet  interval  during 
the  middle  of  the  mission  and  the  various  intensificatimis  semi 
after  Marth  1991. 

Figure  6  quantifies  the  differences  seen  between  the 
measured  and  modeled  flux  profiles.  A  ’difference  ratio*  is 
computed  for  each  L  bin  for  each  day  of  the  missicm.  The 
ratio  is  defined  as  the  absolute  value  of  (measured  flux  - 
modeled  flux)  /  measured  flux.  The  three  curves  in  Figure  6 
are  the  fraction  of  days  (y-axis)  that  the  ’difference  ratio’ 
calculated  from  the  quasi-static  Ap,}  model,  the  NASA 
AE8MAX  model  and  &e  CRKES-static  model,  is  less  than 
unity  for  each  L  (x-axis).  In  gmieral  terms,  the  plots  represent 
the  fraction  of  time  the  models  give  numbers  within  a  factor 
of  2  of  the  measured  values.  Overall,  the  Ap,}  quasi-static 
model  gives  the  best  agreement  and  the  NASA  AE8MAX  the 
worst.  The  Ap,}  models  are  widiin  a  factor  of  2  of  the 
measured  values  ferly  consistently  (SO*  of  the  time)  over  all 
L’s  from  2.5  to  6.5  Rg.  The  NASA  model,  on  the  other 
hand,  shows  serious  problems  for  all  but  the  lowest  L-values. 
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Figufe  6.  Line  bistognm  for  the  ceie  where  the  ’difference  ratio’  between  the  fluxes  of  the  vuiout  models  end  the  measured  data  is  less  than 
1  for  each  L  bin.  The  ’difference  ratio’  is  defined  as  the  absolute  value  of  (measured  flux-modeled  fiux)/meaaured  flux. 


Upcm  exazninatioja  of  individual  spectral  comparisons  (not 
shown),  it  is  serai  that  the  NASA  spectra  are  much  harder  than 
the  CRRES  spectra,  with  the  result  that  NASA  predicts  much 
higher  fluxes  at  higher  energies  up  to  5  MeV  whidi  is 
consistent  with  the  higher  NASA  pre^cted  electron  doses  [8, 
9]. 

In  conclusion,  the  NASA  models  are  the  most  widely  used 
models  for  predicting  electron  radiation  to  near-Earth  space 
systems.  We  have  shown  that  the  AE8MAX  model  has  severe 
deficiraicies  when  compared  to  recent  CRRES  measurements, 
particularly  at  energies  >  2  MeV  and  for  distances  beyond 
3.5  R^.  Trying  to  model  a  highly  dynamic  particle  population 
with  a  single  static  model  will  never  produce  reasonable 
results  under  all  conditions.  As  a  minimum,  two  models, 
representing  quiet  and  active  magnetospheric  conditions,  are 
needed.  The  approach  presraited  here,  driven  by  the  magnetic 
activity  index  Ap,  provides  a  promising  methodology  for 
developing  a  fully  functional  quasi-static  outer  zone  electron 
model.  By  supplemraiting  the  current  CRRES  database  with 
future  data  sets  to  inqnove  model  statistics,  this  quasi-static 
model  could  become  the  standard  for  predicting  electron  fluxes 
(dose)  for  future  sqiace  missions. 
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